Tunnelling couplings in discrete lattices, single particle band structure and eigenstates 

of interacting atom pairs 
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By adjusting the tunnelling couplings over longer than nearest neighbor distances it is possible in 
discrete lattice models to reproduce the properties of the lowest energy band of a real, continuous 
periodic potential. We propose to include such terms in problems with interacting particles and 
we show that they have significant consequences for scattering and bound states of atom pairs in 
periodic potentials. 
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INTRODUCTION 



Spatially periodic potentials induced by off-resonant 
laser fields have been applied in numerous studies to 
modify the dynamics of ultra-cold bosonic gases. In 
the "mean-field" regime with many atoms, the potential 
changes the transport properties and effective mass of the 
atoms [H @, Q j m a dilute gas the lattice enables tuning 
of the ratio between kinetic and interaction energy and 
observation of a Tonks gas [1, [B[ , and in the limit of well 
confined atoms in deep lattice wells collisionally induced 
collapses and revivals of the matter wave mean field has 
been observed [f| . The observation of the transition from 
a superfluid to a Mott insulator is a hallmark experiment 
on a zero temperature phase transition made possible by 
the use of an optical lattice @, H, [s[ ■ 

To deal properly with the many-body aspects of the 
quantum state, it has been useful to approximate the mo- 
tion in the periodic potential by a discrete lattice model, 
specifying for a single atom the atomic wave function 
only on positions corresponding to the potential minima 
of the true potential. The kinetic energy gives rise to 
tunnelling couplings between the sites, whereas the po- 
tential energy is described by local values of the potential 
on each site. If the potential attains the same value on 
each site, the eigenvectors of the Hamiltonian are char- 
acterized by amplitudes that experience a fixed phase 
rotation as one steps through the lattice. These states 
are (quasi-)momentum eigenstates, and their energies de- 
pend on the tunnelling coupling coefficients. The dis- 
crete lattice model with only nearest neighbor tunnelling 
couplings does not reproduce important features in the 
momentum dependence of energy in the lowest energy 
band of a real potential. In this paper we show that it is 
relatively easy to remedy this problem by incorporation 
of "long distance" tunnelling couplings, and we propose 
that such couplings be incorporated in studies of many- 
body dynamics in periodic potentials, since binding and 
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collisional interactions may depend crucially on the dis- 
persion properties of atoms in the periodic environment. 

As a particular example we consider two interacting 
atoms, and in particular the phenomenon of repulsively 
bound atom pairs in optical lattices [T3, HH, 

12, 131. The 

stability of a repulsively bound pair may be understood 
qualitatively from the observation that the repulsive in- 
teraction gives two atoms an energy that lies within the 
energy band gap of free motion in the periodic poten- 
tial, forcing the atoms to stay together in a discrete 
bound state lying above the energy continuum of scat- 
tering states. The lack of dissipative channels, unlike, 
e.g., the phonons in real crystals, makes the observation 
of repulsively bound atom pairs possible in optical lat- 
tices pjj. 

The paper is structured as follows. In Sec. II, we 
analyze the eigenstates of a single particle in a periodic 
potential, and we show that the lowest energy band can 
be obtained also from a discrete lattice model with suit- 
ably selected tunnelling couplings. In Sec. Ill, we turn 
to the description of two atoms with a short range in- 
teraction, and we identify the continuum of scattering 
states in the lowest Bloch band. In Sec. IV we iden- 
tify the discrete states following from the attractive and 
repulsive binding mechanisms, and we show that these 
states have a number of physical properties that depend 
markedly on our more precise description of the band 
structure of the system. Sec. V concludes the paper and 
discusses possible consequences of incorporating its re- 
sults in the treatment of many-body problems, e.g., via 
the Bose-Hubbard Hamiltonian. 



II. SINGLE-PARTICLE DYNAMICS 

We consider an atom with mass m that moves in an 
off-resonant standing wave laser field with wave length 
A. Due to the light induced energy shift, the atom ex- 
periences a potential V(x) = Vq cos(27rx/a) with period 
a = A/2 and potential strength Vq proportional to the 
light intensity and inversely proportional to the optical 
detuning from atomic resonance. According to Bloch's 
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theorem, the Hamiltonian 

is diagonalized by wave functions (f> n k , which are products 
of plane wave factors exp(ifcir) and periodic functions. 
The energies of these states E n (k) describe energy bands 
as the quasi-momcntum k is varied across the first Bril- 
louin zone [— n/a, ir/a]. It is convenient to represent the 
energy in units of the lattice recoil energy En = 2 m\ i • n 
is the band index, and in the following we shall restrict 
our analysis to the lowest band with n = 1 and omit the 
quantum number n from our equations. An alternative 
basis for the quantum states in the lowest band is the 
spatially localized Wannier functions 

Wj (x) = ^=J2 eikaj ^( x ) ( 2 ) 

V k 

describing atoms localized in the j th well in the lattice (to 
ensure constructive interference at the j th well we need 
a phase convention and choose 4>k{x — 0) to be real and 
positive for all k). We have, for convenience, assumed a 
finite number M of lattice periods and periodic boundary 
condition on the entire system. 
By evaluating all matrix elements 

J d = -H hj+d = -(wj\JF\w j+d ,) 

f f h d 2 \ (3) 

= yd* «*(*)*(-— —+ V (x)j Wd (x), 

where we have used the translational invariance of the 
problem, the Hamiltonian can be explicitly rewritten in 
the Wannier basis, 

J?=-J2 Jd(\«>j)(wj+d\ + \wj+d)(v>i\)' ( 4 ) 

So far no approximations were done except for the omis- 
sion of the higher bands, which can be reintroduced in 
the formalism if needed. The Hilbert space is now rep- 
resented by a discrete basis of localized Wannier states, 
and it is only a matter of re- interpretation to read Eq. Q 
as a discrete lattice model of the problem. 

To interpret J d in terms of the band structure, we in- 
sert the expression for the Wannier functions ((2|) into the 
integral for J d and obtain 

k,q v / 

x <j> k (x)e ikad = ~ Yl E(k)e ikad . (5) 

k 

where we observe that the J d amplitudes are the Fourier 
coefficients of the band structure. The band structure is 
an even function of k and therefore 

E(k) = -Jo -2^2 J d cos(kad). (6) 

d>0 
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FIG. 1: Band structure for a lattice with Vo = 0.5Er. Solid 
line: Exact band structure. Dashed line: Apparent band 
structure in a discrete lattice model including only nearest 
neighbor tunnelling (exact only in the large Vo limit). Dash- 
dotted line: Band structure in a discrete lattice model in- 
cluding both Ji and J2. Insert: The importance of higher Jd 
coefficients as function of the lattice strength Vo- The filled 
circles indicate the Fourier coefficients of the free space kinetic 
energy spectrum in the first Brillouin zone. 



Without loss of generality we can set Jo = in the fol- 
lowing. 

The application of lattice models which include only 
nearest neighbor tunnelling couplings only provide us 
with a simple cosine shaped band structure, as depicted 
in Fig. Q] (dashed line). Including more J d terms, the J d - 
values are compared to J\ in the insert of Fig. [TJ gives 
a more accurate energy dispersion as indicated by the 
dash-dotted line in the figure. As the lattice strength is 
ramped up in the limit of strong potentials Vq » Er, 
the importance of the higher terms decreases and only J\ 
differs significantly from zero. For J2 to be less than 1% 
of Ji, V has to exceed 5.327^ Q. 

We observe that in the opposite limit where Vo — > 0, 
J d coincides with the Fourier transform of the free space 
energy dispersion parabola, i.e. J d = —j^(ad)~ 2 (—l) d , 
with absolute values plotted as filled circles in the insert 
of Fig. [TJ 

Fig. [T] clearly shows significant and well known dif- 
ferences between the real band structure and the simple 
cosine dependence found for nearest neighbor tunnelling 
models. These differences are visible if one tries to per- 
turb and excite motion in the system, and as we shall see 
in the following section they also have important conse- 
quences for the states of interacting atoms in the periodic 
potential. 
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III. ATOMIC SCATTERING STATES IN THE 
LATTICE 

We now proceed to the two-body problem and find the 
eigenstates ip(x,y) of a pair of particles moving on the 
lattice with discrete one-dimensional position coordinates 
x and y, and interacting by short range interactions. The 
Hamiltonian of this system is assumed to be of the form 

* = - E Jd ( A - } + A i d) + 4 ) + U6 *,v w 

where the action of the il d th — Laplacians" on a wave func- 
tion are given by A x f(x) — (f(x + ad) + f(x — ad) — 
2f(x)) with step size ad. The sum over the discrete 
Laplacians incorporates the tunnelling terms ([4]) which 
in turn describe both the kinetic energy and the periodic 
potential energy of the particles, i.e. the full dynamics of 
independent particles in the lattice. The interaction be- 
tween particles should be computed as a matrix element 
between products of Wannier states of the real physical 
interaction, and in case of a short range interaction, rep- 
resented by a Dirac delta function gir>5(x — y), one would 
thus find matrix elements in the product state basis, 

Uijki = 9id J dx Wi(x)*Wj(x)*w k (x)wi(x) (8) 

but we will assume in the following that they fall off 
rapidly for Wannier functions situated in different wells, 
and we apply just the single on-site Kronecker delta func- 
tion interaction in ([7]). 

We will first discuss the scattering states, where the 
atoms are asymptotically free. As in normal scatter- 
ing theory, the continuum part of the Hamiltonian |(7J|, 
and the sum of the single particle Hamiltonians of non- 
interacting particles are unitarily equivalent. This im- 
plies that the energy spectrum is given by simple addi- 
tion of contributions from two identical band structure 
calculations as the ones shown in Fig. 1. We can also 
solve the scattering problem explicitly and find the eigen- 
states, and since we shall need this description to ac- 
count for the bound states, we shall proceed and find the 
two-body wave function in terms of the center-of-mass, 
Z = (x + y)/2, and relative coordinate, z — x — y. Note 
that the problem of particles moving in the full cosine po- 
tential does not separate in these coordinates, but in the 
lattice model, the Hamiltonian does not explicitly depend 
on the shape of the lattice potential wells (the interaction 
is built into the "Laplacian" and the Wannier functions 
making up the discrete states), and separation is possi- 
ble. When we apply the Hamiltonian (0 to the wave 
function Ansatz, ip(x,y) — e tKZ ipK(z), for each value of 
K we have to solve the one-body Schrodinger equation 
for the relative position coordinate, 

[ -2 Y, cos(\Kda)J d (a« + 2) + US Z<0 J ^ K (z) 

V d>0 J 

= E4, K {z). (9) 
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FIG. 2: Bound and scattering state spectra for a lattice with 
Vo = 0.5Er and U = 2J\. The shaded bands show the scat- 
tering states with density of states indicated by the shading, 
(a) Calculation including only Ji. (b) Calculation including 
Ji and J2. (c) Calculation including Ja up to d = 6. 

It is readily verified that plane wave functions of the 
form il)® K (z) — e tkz are eigenfunctions of the unperturbed 

Hamiltonian Hq — —2 J2d>o cos(^Kda) Jd(Az +2) with 
energies 

e K (k) = -4 J^Vd cos(\Kad) cos(kad). (10) 

d>0 

As noted above, the spectrum of scattering states, but 
not the scattering wave functions themselves, are inde- 
pendent of the short range scattering potential, and Eq. 
(HHJ) thus provides also the spectrum of interacting atoms 
on the lattice. This spectrum is displayed in Fig. HJa)- 
(c) with inclusion of different numbers of tunnelling co- 
efficients Jd- The different shades of grey in the figure 
show the density of states. For fixed center-of-mass mo- 
mentum K the density of states has a peak around the 
maximal and minimal energy. This is due to the flat- 
ness of the cos(kad) function in Eq. (jTUJ) near its ex- 
tremal values. A complete degeneracy in k can be ob- 
served due to the front factor cos( ^Ka) in Eq. (fit))) , 
when only nearest neighbor coupling is included and 
K = ±7r/a. This degeneracy, however, is lifted by the 
amount Ae ±7r/a = 8E m (-l)"V 2m in Fig. ED>)-(c). 

Another pronounced feature in the spectra of Fig. 
[2] is the extra van Hove singularity in the scattering 
states for values of the center-of-mass momentum close 
to K = ±7r/a. To see where it comes from we note 
that the non-interaction system separate, so we could 
have chosen a product wave function ^ = ip qi (x)ip g2 (y) 
with total energy E(qi) + E(q2), where qj is the quasi- 
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FIG. 3: Band structure of two-particles epc{k) for different 
values of K (solid line) and the band structure of the indi- 
vidual atoms, i.e. + k) and — k) (dotted lines). In 
(a)-(c) the band structure is approximated by a cosine and in 
(d)-(f) higher terms up to Jq is included in the band structure, 
see Eq. ©. 



momentum of the individual atoms and E is the one- 
body band structure ©■ If we infer the center-of-mass 
quasi-momentum K — q\ + q2 and relative motion quasi- 
momentum k = (q\ — q^/S, we can see that this is indeed 
the same energy ex{k) as derived above 



e K (k) 

(cos(f 



-E( qi 



d>0 

Jdcos 

d>0 



f + k) + E (q 2 = f - k) 
fc) ad + cos (-j 

^ cos had. 



- k) ad) 



(11) 



The extra van Hove singularities can be understood 
from Fig. [3J where the band structure for each atom 
is plotted together with the total energy as a function 
of the relative momentum k. The extra van Hove sin- 
gularity around K = ±ir/a in Fig. EKb)-(c) is due to 
the appearance of a local minimum in the two-body en- 
ergy spectrum at k — ixja. This phenomenon sets in 
ir/a) = as shown in Fig. EJe), where 



when {k 



and 



de K {k) _ 







K = 0.4 n/a. We note that — — — u anu — 
cannot be simultaneously fulfilled when the band struc- 
ture is a pure cosine function, i.e., when only nearest 
neighbor tunnelling is included in the model. When 
K = ±7r/a, as in figure the energy minima at k = 
and k — n/a are degenerate and therefore we see only 
two maxima in the density of state in Fig. [2^c). 

In the next section we include the atom-atom interac- 
tion, V(z) = US z ,q, to account for bound states. 



IV. REPULSIVELY AND ATTRACTIVELY 
BOUND STATES 

The bound states of the system can be obtained from 
the Green's function solution of the Schrodinger equation 
([9]). We first solve the problem in the absence of the 
interaction, 



(E-H )G° K (E,z) = 6 zfi 



(12) 



to obtain the unperturbed Green's function. The full 
Green's function in position space can now be obtained 

G K {E,z) = G° K (E,z) 

+ J dz'G° K (E,z- z')V{z')G K {E,z') 

= G K (E,z) + G° K (E,z)UG K (E,z = Q) 
= G° K (E,z) + G0 ^ Z ^ E > Z = V 



1 - UG° K (E,z = 0) 



G° K (E,z) 



1 - UG° K (E,z = 0)' 



(13) 



Whenever there is a bound state Gk has a pole, there- 
fore the bound state energies are determined solely by 
the interaction strength U and the unperturbed Green's 
function 



l = UG° K (E b ,z = 0). 



(14) 



The unperturbed Green's function, in turn, is read- 
ily obtained from the momentum space representation, 
G° K (E,z) = (27:)- 1 Jdk G° K (E,ky kz , 



G° K (E,k) 



E — e K (k) + irj 



(r?^0+). 



(15) 



The integral of Eq. (|15p can be found analytically if 
e^-(fc) is a simple cosine, and in the general case it can be 
computed numerically and the solution to Eq. (|14[) can 
be found by a simple numerical search. 

We can combine Eqs. (fl~2j) and (fT4| to obtain the re- 
lation 



(E-H )G° K (E b ,z) = 



1 



G° K (E b ,z = 0) 
U8 z , G° K (E b ,z), 



S zfi G° K (E b ,z) 



(16) 



This shows that G° K (E b ,z) is itself a solution to the 
Schrodinger equation, i.e., it provides directly the bound 
state wave function -0^(z), and G° K {E b ,k) provides the 
wave function Tpx(k) in momentum space. 

There are discrete bound states for both positive and 
negative U corresponding to repulsive and attractive in- 
teraction, respectively. The energy of the repulsively 
and attractively bound atom pairs lie above and below 
the scattering continuum as plotted in Fig. [21 The re- 
pulsively bound pair state was recently experimentally 
demonstrated by Winkler et al. [Tol | . 
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FIG. 4: Repulsively bound states obtained with nearest neigh- 
bor interaction only. First column shows i/)K{k) and the sec- 
ond tpK(z) both obtained from Eqs. (Tl5l fT6|l . see text for 
details. 
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FIG. 5: Repulsively bound states obtained with six nearest 
neighbor interactions included. First column shows ipK(k) 
and the second ipii(z) both obtained from Eqs. (|15l I16p . see 
text for details. 



If we only include J\ we have the following symmetry 
in the continuum spectrum, Eq. (|10|) . 

e K (k) = -e K (k + Tr/a). (17) 

This implies, according to Eq. (fT5|) . when E is outside 
the scattering continuum, 

G° K (E,k) = -G° K (-E,k + n/a), (18) 

and hence 

1 =UG° K (E ab (K),z = 0) 
=(-U)G° K (-E ab {K),z = 0) (19) 
=(-U)G° K (E rb (K),z = 0), 

which relates the repulsively (rb) and attractively (ab) 
bound atom pair energies: 

E ab (K) = -E rb (K). (20) 

Eq. P^|) now implies, if we include only Ji, that the 
attractively and repulsively bound wave functions are re- 
lated by 

^(k)=^(k + 7r/a). (21) 

This symmetry is also pointed out in 12], but it is 
important to note that when the higher Jd tunnelling 
couplings are included, Eq. (|17|) and hence the following 
equations ([2TJI l2"Tj) are no longer valid. 

The shape of the band structure is altered when Jd>i 
are included, but more importantly the wave functions 
are altered significantly, compare Figs. 0] and 



In Fig. [4] only Jd=i is included. As one would expect 
the bound states are composed of the quasi-momentum 
components lying closest in energy to the bound state. In 
the region around K = this correspond to relative mo- 
mentum k = ir/a and the formation of bound pairs can 
in this region be explained by the fact that the peak in 
fc-distribution corresponds to the individual atoms hav- 
ing quasi-momentum q\ = qi = ir/a. Here the effective 
mass is ne gati ve and therefore a repulsive force causes 
attraction [13 |. As K approaches ±ir/a the scattering 
states with different relative momenta become degenerate 
in energy and the momentum distribution of the bound 
state goes from a peaked function Fig. [Ha) to a constant 
function Sfe) , which corresponds to a delta function in 
relative coordinate, H{f). 

In Fig. O Jd up to d < 6 has been included in obtaining 
the fc- and r-distributions of the repulsively bound states. 
As discussed previously; when the extra van Hove singu- 
larity sets in Fig. [UJe), the two-body spectrum broadens 
around k — n/a and this can also be observed in the de- 
composition of the bound states fig. EJc). In Fig. [SJe) we 
clearly see the consequences of lifting the fc-degeneracy 
at K — ztn/a. The momentum distribution is no longer 
constant but instead peaked around the two momenta 
k = ±7r/2a which compose the top of the scattering band, 
as can be seen in Fig. [SJf). 

In the experiment by Winkler et al. the atoms 

were prepared in the K = bound state. We suggest that 
by exposing the system to a force field the bound pairs 
can be prepared in any instate. This can be realised 
in practice by chirping the lattice lasers so the lattice 
potential is accelerated in the laboratory frame [l5| ■ To 
investigate the effect on the two-atom states, we make 
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FIG. 6: Illustration of the tunnelling and interaction coeffi- 
cients in the Bose-Hubbard Hamiltonian with "long distance" 
tunneling. 

the new time dependent Ansatz 

*(x,y,t) = e iK ^i>(x-y,t) (22) 

with a time dependent center-of-mass quasi-momentum, 
K{t). Inserting this Ansatz into the time dependent 
Schrodinger equation with the Hamiltonian Q and the 
linear term due to the force, — F ■ (x + y), we obtain 

.d., . 
in—ip{z,t) = 

(-2^2 cos{\K{t)da)J d (AM + 2) + US Z . \ i/,(z, t) 

(23) 

if we choose 

IF 

K(t) = — . (24) 

Since there is an energy gap between the bound state 
and the scattering states, we may apply the adiabatic 
approximation to Eq. (|23p and simply replace ij}{x — 
y,t) by the eigenstate iPk{x — y) of Eq. © with K = 
K{t). Provided we accelerate the lattice slowly enough 
we can thus explore all the wave functions ipK (x — y) of 
the system with K — 2Ft/h. 

V. DISCUSSION 



in this paper shown that the lowest band in the single 
particle band structure in such a potential can be de- 
scribed very accurately with a discrete lattice model with 
tunnelling couplings to a sufficiently high number of re- 
mote sites. Such modelling may be particularly impor- 
tant when one studies scattering processes where energy 
and momentum are conserved. We have previously stud- 
ied the break-up of a single condensate as a four wave 
mixing process [16[, which cannot take place if the band 
structure is a simple cosine function. In the present work 
we described the continuum of scattering states of atoms 
interacting by a short range attractive or repulsive inter- 
action, and we identified the two-atom repulsively and 
attractively bound states. 

Simple models for the periodic potentials may be of 
particular interest when systems with a large number 
of particles are studied. In this case, one has recourse 
to second quantization, and the restriction to a discrete 
set of Wannier states is a first step towards an accu- 
rate treatment of the many-body properties, e.g., how 
the atoms distribute themselves among potential wells. 
See ref. for a whole list of projects for further 

studies. So far many problems were dealt with by the 
Bose-Hubbard Hamiltonian, @, Q [13, H 0, H3, ei- 
ther solved exactly numerically or exposed to further 
simplifying approximations. Based on the experience of 
the present work, we suggest that in cases where col- 
lisions with momentum and energy conservation play a 
role, one should extend the conventional Bose-Hubbard 
Hamiltonian with only nearest neighbor tunnelling to in- 
clude also tunnelling terms over longer distances, Jf = 
~ ^iM>o Jd(&]aj + d+a^ +d &j)+Udjajajaj, as depicted in 
Fig. [6j This modified Hamiltonian, or further extended 
with non-local interaction terms [2l|, [22j , would seem to 
be a good starting point for detailed studies of collisions 
and instabilities in moving condensates in periodic poten- 
tials. We are currently pursuing such studies, including 
also the interplay of repulsive binding and real molecular 
binding of atoms in optical lattices. 



The states and energies of particles in periodic poten- 
tials depend on the details of these potentials. We have 
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